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Abstract 

We give a survey of the two remarkable analytical problems of 
quantum information theory. The main part is a detailed report of the 
recent (partial) solution of the quantum Gaussian optimizers problem 
which establishes an optimal property of Glauber’s coherent states - 
a particular instance of pure quantum Gaussian states. We elaborate 
on the notion of quantum Gaussian channel as a noncommutative 
generalization of Gaussian kernel to show that the coherent states, 
and under certain conditions only they, minimize a broad class of 
the concave functionals of the output of a Gaussian channel. Thus, 
the output states corresponding to the Gaussian input are “the least 
chaotic”, majorizing all the other outputs. The solution, however, is 
essentially restricted to the gauge-invariant case where a distinguished 
complex structure plays a special role. 

We also comment on the related famous additivity conjecture, 
which was solved in principle in the negative some five years ago. 
This refers to the additivity or multiplicativity (with respect to tensor 
products of channels) of information quantities related to the classical 
capacity of quantum channel, such as (1 —)• p)-norms or the minimal 
von Neumann or Renyi output entropies. A remarkable corollary of 
the present solution of the quantum Gaussian optimizers problem is 
that these additivity properties, while not valid in general, do hold in 
the important and interesting class of the gauge-covariant Gaussian 
channels. 
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1 Introduction 

The quantum Gaussian optimizers problem is an analytical problem that 
arose in quantum information theory at the end of past century, and which 
has an independent mathematical interest. Only recently a solution was 
found [23], [33] in a considerably common situation, while in full generality 
the problem still remains open. To explain the nature and the difficulty of 
the problem we start from the related classical problem of Gaussian maximiz¬ 
ers which has been studied rather exhaustively, see Lieb [50] and references 
therein. Consider an integral operator G from Lp (M*) to Lg (M^) given by a 
Gaussian kernels (i.e. exponential of a quadratic form) with the {q ^ p) — 
norm 

IIGIIg^p = sup IIG/II /ll/ll = sup \\Gf\\ (1) 

m ll/ll <1 
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Under certain broad enough assumptions concerning the quadratic form 
dehning the kernel, and also p and g, this operator is correctly dehned, and 
the supremum in ([1]) is attained on Gaussian /. Moreover, under some addi¬ 
tional restrictions any maximizer is Gaussian. As it is put in the title of the 
paper [50]: “ Gaussian kernels have only Gaussian maximizers”. 

Knowledge that the maximizer is Gaussian can be used to compute exact 
value of the norm ([T]); in fact a starting point of the classical Gaussian max¬ 
imizers works were the result of K.I. Babenko [5] and a subsequent paper 
of Beckner [6] which established the best constant in the Hausdorff-Young 
inequality concerning the {p ^ p') —norm, {p~^ + 1 < P < 2), of 

the Fourier transform (which is apparently given by a degenerate imaginary 
Gaussian kernel). 

A difficulty in the optimization problem ([1]) is that it requires maximiza¬ 
tion of a convex function, so the general theory of convex optimization is not 
of great use here (it only implies that a maximizer of ||G/||p belongs to a face 
of the convex set \\f\\q < 1). Instead, the solution is based on substantial 
use of the classical Minkovski’s inequality and the related multipicativity of 
the classical {q — )■ p) —norms with respect to tensor products of the integral 
operators. 

A notable application of these classical results to a problem in quantum 
mathematical physics was Lieb’s solution [51] of Wehrl’s conjecture [63]. Let 
p be a density operator in a separable Hilbert space "H representing state of 
a quantum system; the “classical entropy” of the state p is dehned as 0 

f z 

Hci{p) = - / Pp{z) logpp{z) -, 

Jc ^ 

where Pp{z) = {z\p\z) is the diagonal value of the kernel of p in the system 
of Glauber’s coherent vector^ ^ [S], [33]- The conjecture was 

that Hci{p) has the minimal value if p is itself a coherent state i.e. projec¬ 
tor onto one of the coherent vectors. Lieb [51] used exact constants in the 
Hausdorff-Young inequality for Lp-norms of Fourier transform 151.0 and the 
Young inequality for convolution [6] to prove similar maximizer conjecture 
for f[x) = x'P and considered the limit limp^i(l — p)~^ (1 — x^) = — x logx. 

^ Throughout the paper the base of logarithm is a fixed number a > 1. In information 
theory the natural choice is a = 2, then all the entropic quantities are measured in “bits”. 

^ In analysis, they correspond to complex-parametrized Gaussian wavelets. Notice 
that this is the only place in the present article where we formally used Dirac’s notations, 
uncommon among mathematicians. 
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Recently, Lieb and Solovej [S2], by using a completely different approach 
based on study of the spin coherent states, strengthened the result of [SI] 
by showing that the coherent states minimize any functional of the form 
/c/(Pp(^))^! where f{x),x G [0,1] is a nonnegative concave function with 

/(O) = 0. 

In the language of quantum information theory, the affine map G : p ^ 
Pp{z), taking density operators p (quantum states) into probability densities 
Pp{z) (classical states), is a “ quantum-classical channel” [39]. Moreover, it 
transforms Gaussian density operators p (in the sense dehned below in Sec. 
O) into Gaussian probability densities, and in this sense it is a “ Gaussian 
channel” . From this point of view, Wehrl entropy Hci{p) is the output 
entropy of the channel, and Lieb’s result says that it is minimized by pure 
Gaussian states p. Moreover, the corresponding result for f{x) = x^ can be 
interpreted as “ Gaussian maximizer” statement for the norm ||G||^^p . Notice 
that the case q = 1, which is excluded in the classical problem for obvious 
reasons, appears and is the most relevant in the quantum (noncommutative) 
case. 

The quantum Gaussian optimizers problem described in the present paper 
refers to Bosonic Gaussian channels - a noncommutative analog of Gaussian 
Markov kernels and, similarly, requires maximization of convex functions (or 
minimization of concave functions, such as entropy) of the output state of 
the channel, while the argument is the input state. A general conjecture 
is that the optimizers belongs to the class of pure Gaussian states. The 
conjecture, hrst formulated in [12] in the context of quantum information 
theory, however natural it looks, resisted numerous attacks for several years. 
Among others, notable achievements were the exact solution for the classical 
capacity of pure loss channel [21] and a proof of additivity of the Renyi en¬ 
tropies of integer orders p [21] for special channels models. Even restricted 
to the class of Gaussian input states, the optimization problem turns out to 
be nontrivial [56], [31]. There was some hope that in solving the problem, 
similarly to Wehrl’s conjecture, one could also use the classical “ Gaussian 
maximizers” results. However the solution found recently by Giovannetti, 
Holevo, Garcia-Patron [23], and Mari, Giovannetti, Holevo [53] uses com¬ 
pletely different ideas based on a thorough study of structural properties of 
quantum Gaussian channels. As it was mentioned, a solution of the classi¬ 
cal problem uses the Minkowski inequality and the implied multiplicativity 
of {q —)■ p)-norms. However, the noncommutative analog of the Minkowski’s 
inequality [T2] is not powerful enough to guarantee the multipicativity of 
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norms (or additivity of the corresponding entropic quantities). Moreover, 
the related long-standing additivity problem in quantum information theory 
[3^ was recently shown to have negative solution in general [26]. We show 
that, remarkably, a solution of the quantum Gaussian optimizers problem 
given in [23] implies also a proof of the multipicativity/additivity property 
in the restricted class of gauge-covariant or contravariant quantum Gaussian 
channels. 

It would then be interesting to investigate a possible development of 
such an approach to obtain noncommutative generalizations of the classi¬ 
cal “ Gaussian maximizers” results for {q —)• p) —norms. Such generalization 
could shed a new light to the hypercontractivity problem for quantum dynam¬ 
ical semigroups and related noncommutative analogs of logarithmic Sobolev 
inequalities, see e.g. [62] . 

2 The additivity problem for quantum chan¬ 
nels 

2.1 Definition of channel 

Let "H be a separable complex Hilbert space, -C('H) the algebra of all bounded 
operators in H and T('H) the ideal of trace-class operators. The space T('H) 
equipped with the trace norm H-Hj^ is Banach space, which is useful to consider 
as a noncommutative analog of the space Li. The convex subset of T('H) 

&{n) = {p : p* = p > 0, Tip = 1} , 

is a base of the positive cone in Operators p from G{'H) are called 

density operators or quantum states. The state space is a convex set with 
the extreme boundary 

^i(n) = {p : p > 0, Trp = 1, p^ = p} . 

Thus extreme points of which are called pure states, are one-dimensional 

projectors, p = P,p for a vector -0 G "H with unit norm, see, e.g. [55]. 

The class of maps we will be interested is a noncommutative analog of 
Markov maps (linear, positive, normalized maps) in classical analysis and 
probability. Let be the two Hilbert spaces, which will be called input 

and output space, correspondingly. A map $ : T('H^) —)■ T(?fs) is positive 
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ifX > 0 implies > 0, and it is completely positive j61j . j54j if the maps 
$ (g) Id(rf) are positive for all d = 1, 2,, where Id(d) is the identity map of 
the algebra 2,d = of complex d x d— matrices. Equivalently, for every 

nonnegative dehnite block matrix ^ the matrix ^ is 

nonnegative dehnite. 

A linear map <h is trace-preserving if Tr<f)[X] = TrX for all X G %{'Ha)- 

Definition Quantum channel is a linear completely positive trace-preserving 
map <h : %{1-La) —t Letter A will be always associated with the input 

of the channel, while B with the output. Sometimes, to abbreviate notations, 
we will write simply $ : A —)■ S. ■ 

Apparently, every channel is a positive map taking states into states: 
<h[©('Hyi)] C Since T('H) is a base-normed space, this implies [IT] 

that $ is a bounded map from the Banach space %{'Ha) to The dual 

of the map <h is uniquely dehned by the relation 

Tr$[X]F = TrX<h*[F]; X e ^IUa), Y e ^{Hb), (2) 

and it is called dual channel. The dual channel is linear completely positive 
*—weakly continuous map from 2,{l-iB) to SYQdA), which is unital. = 

Here and in what follows I with possible index denotes the unit operator 
in the corresponding Hilbert space. 

There are positive maps that are not completely positive, a basic example 
provided by matrix transposition X —)• in a hxed basis. 

From the dehnition of complete positivity one easily derives [39] that 
composition of channels $2 ° 'Ll dehned as 

^2 0^l[X] = <h2[<hi[X]], 

and naturally dehned tensor product of channels 

$1 (g) $2 = ($1 0 Id2) o (Idi 0 $ 2 ) 


are again channels. 


2.2 Stinespring-type representation 

The notion of completely positive map was introduced by Stinespring [61] in 
a much wider context of C*-algebras. This allows also to cover the notion of 
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hybrid channel where the inpnt is qnantnm while the ontpnt is classical or 
vice versa. An example of snch channel was mentioned in Sec. [U We will 
not pnrsne this topic fnrther here, see [39], bnt only mention that complete 
positivity rednces to positivity in snch cases. 

Motivated by the famons Naimark’s dilation theorem, Stinespring estab¬ 
lished a representation for completely positive maps of C*-algebras which in 
the case of qnantnm channel rednces [59] to 

Proposition 1 Let $ : A —)■ i? he a channel. There exist a Hilbert space 
He and an isometric operator V : Ha Hb ®He, such that 

<^[p] = TteVpV*- pe%{HA), (3) 

where Tie denotes partial trace with respect to He- The representation ^ 
is not unique, however any two representations with Vi : Ha Hb ® He^ 
and V 2 : Ha Hb ® He2 are related via partial isometry W : Hei He2 
such that V 2 = {Ib ® hh) hi and Vi = {Ib ® W*) V 2 . 

Consider a representation ([3D for the channel $; the complementary chan¬ 
nel |37|, [H] is then dehned by the relation 

^p]=TibVpV*- peT:{HA). (4) 

From the relation between the different representations (|3D, it follows that the 
complementary channel is nniqne in the following sense: any two channels 
<Fi,<F 2 complementary to $ are isometrically eqnivalent in the sense that 
there is a partial isometry W : He^ ^ He 2 snch that 

^2[p] = W^^[p]W*, <li[p] = W*^2[p\W, (5) 

for all p. It follows that the initial projector W*W satishes $i[p] = W*W^i[p\, 
i.e. its snpport contains the snpport of <Fi[p], while the hnal projector WW* 
has similar property with respect to *h 2 [p]- The complementary to comple¬ 
mentary can be shown isometrically eqnivalent to the initial channel, so that 
<F, $ are called mntnally complementary channels. 

In general, we will say that two density operators p and a (possibly acting 
in different Hilbert spaces) are isometrically equivalent if there is a partial 
isometry W snch that p = WaW*, a = W*pW. Apparently, this is the case 
if and only if nonzero spectra (connting mnltiplicity) of the density operators 
p and a coincide. We denote this fact with the notation p ^ a. We have jnst 
shown that *hi[p] ~ ^ 2 [p] for arbitrary p. 


7 


Lemma 2 Let ^ be a complementary channel ^), then $[P^] ~ $[P^] for 
all e Ha- 

Proof. Let V : 'Ha —)■ Pb ® Pe be the isometry from the representations 
([3]), dl]), then pbe = VPpV* is a pure state in Hb ® He, and the statement 
follows from a basic result in quantum information theory (“Schmidt decom¬ 
position”): if Pbe is a pure state in Hb®He and pb = Tr^p^^, pe = PtbPbe 
are its partial states, then pB ~ pE (see e.g. Proposition 3 in |31]) ■ 

A different name for channel is dynamical map - in nonequilibrium quan¬ 
tum statistical mechanics they arise as irreversible evolutions of an open 
quantum system interacting with an environment [39]. Assume that there is 
a composite quantum system AD = BE m. the Hilbert space 

P = Ha ® He — Hb ® He, (6) 

which is initially prepared in the state pa ® Pd and then evolves according 
to the unitary operator U. Then the output state pb depending on the input 
state Pa = P is 

^b[p] = T^teU{p®pd)U\ ( 7 ) 

while the output state of the “environment” E is the output of the channel 

^e[p\ =P"^bU{p ® pd)U* . ( 8 ) 

If the initial state of D is pure, po = Ppm then by introducing the isometry 
V : Ha Hb ®Hd, which acts as 

Vf) = U{'4) ®'4 )d), ipeHA, 

we see that the relations ([7j), ([8]) convert into ([3]), (jl]), and ^e is just the com¬ 
plementary of Notice also that both partial trace and unitary evolution 
are completely positive operators, hence the maps (|7]), (|HD are completely 
positive; vice versa, any quantum channel has a representation of such a 
form, see, e.g. [39] . 

Vast literature is devoted to study of quantum dynamical semigroups 
(noncommutative analog of Markov semigroups) and quantum Markov pro¬ 
cesses. Stinespring-type representation Q underlies dilations of quantum 
dynamical semigroups to the unitary dynamics of open quantum system in¬ 
teracting with an environment [HI. ra- 


2.3 Entropic quantities and additivity 

Consider the norm of the map $ dehned similarly to ([T]); 

||d>IK^^ = snp||$[X]|m|X||,= snp ||<h[X]||^, (9) 

||X|P<1 

where H-Hp is the Schatten p—norm [55]. As shown in |1], 

||$|ILp= snp Tr<h[p]?’= snp Tr<h[P^]^ (10) 

P&G{Ha) i’GHA 

where the second eqnality follows from convexity of the fnnction x^,p > 1.. 

The qnantnm Renyi entropy of order p > 1 of a density operator p is 
dehned as 

^p(p) = y^logTr/=-^log||p|| (11) 

1 — p 1 — p ^ 

Dehne the minimal output Renyi entropy of the channel <h 

«.(■*>)= taf = (12) 

P&&(n) i — p 

and the minimal output von Neumann entropy 

mt ( 13 ) 

p&eiv.) 

In the limit p ^ 1 the qnantnm Renyi entropies monotonely nondecreasing 
converge to the von Nenmann entropy 

lim Rp{p) = —Trplogp = H{p). 

p—>-i 

In hnite dimensions the set of qnantnm states is compact, hence by Dini’s 
Lemma the minimal ontpnt Renyi entropies converge to the minimal ontpnt 
von Nenmann entropj|§. 

Mnltiplicativity of the norm ([9|) for some channels 'Ll, $ 2 ; 

||$l«)$2|li^p=||<hl|ll^p-||®2|li^p (14) 

is eqnivalent to the additivity of the minimal ontpnt Renyi entropies 

Pp(^>i(g)$2) = Rp(4>i) + Rp(4>2). (15) 

^The corresponding statement is not valid for infinite-dimensional channels (even for 
classical channels with countable set of states), M. E. Shirokov, private communication. 
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Closely related is the similar property for the minimal output von Neumann 
entropy: 

=i^($i) + ^(<h2). (16) 

In hnite dimensions, the validity of flT^ for certain channels ‘hi ,$2 and p 
close to 1 implies (IThD for these channels. 

In the last two relations the inequality < (similarly to the inequality > 
in fll4p l is obvious because the right-hand side is equal to the inhmum over 
the subset of product states p = pi® p 2 - On the other hand, existence of 
“entangled” pure states which are not reducible to product states, is the 
cause for possible violation of the equality for quantum channels. 


2.4 The channel capacity 

The practical importance of the additivity property (fT6|) is revealed in con¬ 
nection with the notion of the channel capacity. To explain it we assume 
that 'Hai'Hb are hnite dimensional for the moment. 

For a quantum channel $, a noncommutative analog of the Shannon 
capacity, which we call x~capacity, is dehned by 




sup 


H $ 






(17) 


where the supremum is over all quantum ensembles, that is hnite collections of 
states {pi ,..., pn} with corresponding probabilities {tti, ..., tt^}. The quan¬ 
tity flTTP is closely related to the capacity C(<I>) of quantum channel <I> for 
transmitting classical information [3l]. The classical capacity of a quantum 
channel is dehned as the maximal transmission rate per use of the channel, 
with coding and decoding chosen for increasing number n of independent 
uses of the channel 

= $ (g) • • • (g) $ 

n 

such that the error probability goes to zero as n —)■ cxo (for a precise dehnition 
see [39]). A basic result of quantum information theory, HSW Theorem [3^ . 
says that such dehned capacity C(<h) is related to C'^(*h) by the formula 

C($) = lim 

n^oo 


10 




Since C^($) is easily seen to be super additive (i. e., + 

C^($ 2 ) ), one has However if the additivity 

<8) $ 2 ) = CJ-hi) + ^^$ 2 ) (18) 

holds for a given channel $1 = $ and an arbitrary channel $ 2 , then 

C,(4®") = nC,(4>), (19) 

implying 

C($)=C,($). (20) 

The reason for possible violation of the equality here, as well as in the cases 
dn, disi), m, is existence of entangled states, which are not reducible to 
product states, at the input of tensor product channel 

2.5 Main conclusions 

Thus it was natural to ask: does the the additivity property ffTOj) holds 
globally, i.e. for tensor product of any pair of quantum channels $ 1 , $ 2 ? 
The problem can be traced back to [S], see also [M]. Quite remarkably, 
Shor [60], see also [19], had shown the equivalence of the global properties of 
additivity of the y— capacity and of the minimal output entropy. 

Theorem 3 J6^ The properties and IfTU) are globally equivalent in the 
sense that if one of them holds for all channels $ 1 , $ 2 , then another is also 
true for all channels. 

The additivity is proved rather simply for all classical channels (see e.g. 
[T5]). but in the quantum case the question remained open for a dozen of 
years, and was ultimately solved in the negative. 

The detailed history of the problem up to 2006 can be found in [M], and 
here we only sketch the basic steps and the hnal resolution. In [T] it was 
suggested to approach the additivity property flTB]l via multiplicativity flTT|) 
of the (1 ^ p) —norms (equivalent to additivity (IT^ of the minimal output 
Renyi entropies). The first explicit example where this property breaks for 
d = dim'H > 3 and large enough p was transpose-depolarizing channel [6^ : 

<i.(p) = ^[/Ti'p-y], (21) 
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where p G is a matrix and its transpose. In particnlar, (IT^ with 

$1 = $2 = fails to hold for p > 4, 7823 if d = 3 (nevertheless, the addi¬ 
tivity of i7(*h) and of Cx(‘^) holds for this channel). Five years later came 
important hndings of Winter |65] and Hayden [28], see also [29], who showed 
existence of a pair of channels breaking the additivity of the minimal ontpnt 
Renyi entropy for all valnes of the parameter p > 1. The method of these 
and snbseqnent works is random choice of the channels, which for hxed di¬ 
mensions are parametrized by isometries V in the representation ([3|), as well 
of the inpnt states of the channels, combined with snfficiently precise proba¬ 
bilistic estimates for the norms (fTOj) . For finite dimensions the corresponding 
parametric sets are compact, and one usually takes the uniform distribution. 
Basing on this progress, Hastings [25] gave a proof of existence of channels 
breaking the additivity conjecture ffT6|) corresponding to p = 1, in very high 
dimensions. Moreover, the probability of violation of the additivity tends to 
1 as the dimensionalities tend to inhnity. Hastings gave only a sketch, and 
the detailed proof following his approach was given by Fukuda, King and 
Moser [18], and further simplihed by Brandao and M. Hordecki uni. Later 
Szarek et ah [5] proposed a proof related to the Dvoretzky-Mirman theorem 
on almost Euclidean sections of high-dimensional convex bodies. 

Although, combined with theorem [3] this gives a definite negative answer 
to the additivity conjectures, several important issues remaine open. All the 
proofs use the technique of random unitary channels or random states and as 
such are not constructive: they prove only existence of counterexamples but 
do not allow to actually produce them. Attempts to give estimates for the 
dimensions in which nonadditivity can happen based on Hastings’ approach 
has led to overwhelmingly high values: the detailed estimates made in [T8] 
gave d ~ 10^^ breaking the additivity by a quantity of the order 10“^. The 
best result in this direction obtained in [7] states that “violations of the 
additivity of the minimal output entropy, using random unitary channels 
and a maximally entangled state state, can occur if and only if the output 
space has dimension at least 183. Almost surely, the defect of additivity is 
less than log 2, and it can be made as close as desired to log 2”. 

While this does not exclude possibility of better estimates, based per¬ 
haps on a different (but yet unknown) models, it casts doubt onto Ending 
concrete counterexamples by computer simulation of random channels. From 
this point of view, the following explicit example given in [25] is of interest. 


12 


Consider the completely positive map 

p d).[p] = Tt 2 P-pP -, pe%{P®H), 

where P_ is the projector onto the antisymmetric subspace 'H_ of which 

has the dimensionality , and the partial trace is taken with respect to 
the second copy of "H. Its restriction to the operators with support in the 
subspace 'H_ is trace preserving, hence it is a channel. It can be shown [39] 
that <h_ = -where I** is the dual to the complementary of the channel 

fl2T]) . For this simple channel the minimal Renyi entropies are nonadditive 
for all p > 2 and sufficiently large d, but unfortunately it is not clear if it 
could be extended to the most interesting range p > 1. 

Coming back to arbitrary channels, it remains unclear what happens 
in small dimensions: perhaps the additivity still holds generically for some 
unknown reason, or its violation is so tiny that it cannot be revealed by 
numerical simulations. This is indeed surprising in view of the fact that the 
physical reason for nonadditivity is entanglement between the inputs of the 
parallel quantum channels, see [39] for more detail. 

On the other hand, these results stress the importance of continuing ef¬ 
forts to hnd special cases where the additivity holds for some reason, and 
can be established analytically. 

A survey of the main classes of such “ additive” channels acting in hnite 
dimensions was presented in [34] ; below we briefly list the most important 
classes of channels d* for which the additivity properties (HE]), (ITSll and (|T5il 
for p > 1 were established with $ = <I)i and arbitrary $ 2 . 

• Qubit unital channels, i.e channels <h : £2 -^2 satisfying <F[J] = I 

[46] . Strikingly, there is still no analytical proof of the additivity for 
nonunital qubit channels, in spite of a convincing numerical evidence 

ra- 

• Depolarizing channel in 2,^ : 

I (P 

4>[p] = (1-p)p + p-Trp, 0 < p < ^ , 

which is the only unitarily-covariant channel, and can be regarded as 
noncommutative analog of completely symmetric channel in classical 
information theory [15]. The additivity properties flT6|l . flT5|) . flT8|l were 
proved by King BU¬ 


IS 





• Entanglement-breaking channels. In finite dimensions these are chan¬ 
nels of the form 

3 

where {Ma} is a resolution of the identity in "Ha: Ma > 0 , X] = 

3 

Ia, and pb G G{Hb) (see [13]). 

For the finite-dimensional entanglement-breaking channels the additiv¬ 
ity of the minimal output von Neumann entropy and of the x~capacity 
was established by Shor [59] and the additivity of the minimal out¬ 
put R enyi entropies - by King [15|. The additivity properties of 
entanglement-breaking channels were generalized to infinite dimensions 
by Shirokov [SB]- 

• Complementary channels. 

The additivity of the minimal output entropy is equivalent for a chan¬ 
nel <I> and its complementary d), see Lemma HI below. The class of 
channels complementary to entanglement-breaking contains the Schur- 
multiplication maps of matrices p = [cjk] j,k=i,...,d in S.d'- 

where [ 7 ^*;] is a nonnegative definite matrix such that 7 ^^ = 1 . 

For these channels, which are also called “Hadamard channels” the 
additivity of the y—capacity was also established [IB]. 

In the next Sections we consider Bosonic Gaussian channels which act in 
infinite-dimensional spaces. One of the main goals of the present paper is to 
show that the additivity holds for a wide class of gauge co- or contravariant 
Gaussian channels, i.e. those which respect a fixed complex structure in the 
underlying symplectic space. 

2.6 Majorization for quantum states 

From now on we again allow the Hilbert spaces in question to be infinite¬ 
dimensional. Denote by ^ the class of real concave functions / on [0,1], such 
that /(O) = 0. For any f E ^ and for any density operator p we can consider 
the quantity 

W(P) = /(A,), 

3 
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where Xj are the (nonzero) eigenvalues of the density operator p, counting 
multiplicity. Note that this quantity is defined unambiguously with values in 
(— 00 , 00 ]. This follows from the fact that f{x) > cx, where c = /(I), hence 
Tr/(p) > cTrp = c. We also will use the fact that the functional p —)■ Tr/(p) 
is (strictly) concave on &{%) if / is (strictly) concave (see e.g. [TT]b 

Denote by Aj(p) the eigenvalues of a density operator p,counting multi¬ 
plicity, arranged in the nonincreasing order. One says that density operator 
p majorizes density operator a if 

k k 

5^Aj(p) > 5^Aj(a), k = l,2,... 

i=i j=i 

A consequence of a well known result, see e.g. m, is that this is the case if 
and only if Tr/(p) < Tr/(cr) for all / G 5^. 

For a quantum channel $ we introduce the quantity 

where the second equality follows from the concavity of the functional p —>■ 
Tr/($[p]) on Moreover, for strictly concave /, any minimizer is of the 

form Pp for some vector ^ G "H. 

In particular, taking f{x) = —x\ogx and f{x) = —we obtain /(<F) = 
H{^) and/(«F) = -|l«F||^^^. 

Lemma 4 For complementary channels, /($) = /($). Hence 114*11^^^ = 
||4)||i^p, H{^) = H{^), Rp{^) = Rp{^), and the multiplicativity (\T^, as 
well as the additivity of the minimal output entropies holds simul¬ 

taneously for pairs of channels 4)i, 4>2 and $ 1 , 4>2. 

Proof. From Lemma El ^[Pp] and 4>[P^] have identical nonzero spectrum 
mPp] ~ <F[Pp]) . Then 

TV/('S>l^’*l) = TrfmP^]) (23) 

since /(O) = 0. Using second equality in fl22|) implies /(4>) = /($). 

The statement about multiplicativity (additivity) then follows from the 
fact that the channel 4)i (g) $2 is complementary to 4>i (g) 4>2. ■ 
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Quantum Gaussian systems 

3.1 Gaussian states and channels 


A real vector space Z equipped with a nondegenerate skew-symmetric form 
A(z,z') is called symplectic space. In what follows Z is hnite-dimensional, 
in which case its dimensionality is necessarily even, dimZ = 2s [IH]- A basis 
{ej, hj]j = 1 ,..., s} in which the form A{z, z') has the matrix 


A = diag 


0 

-1 



(24) 


is called symplectic. The Weyl system in a Hilbert space "H is a strongly 
continuous family {W{z);z G Z} of unitary operators satisfying the Weyl- 
Segal canonical commutation relation (CCR) 


W{z)W{z') = exp[—^A{z, z')]W{z + z'). (25) 

Thus z ^ W{z) is a. projective representation of the additive group of Z. 
We always assume that the representation is irreducible. The Stone-von 
Neumann uniqueness theorem says that such a representation is unique up 
to unitary equivalence. It is well-known, see e.g. [55], that there is a family 
of selfadjoint operators z -R(^) with a common essential domain V such 


that 


W{z) = expiR{z), 


moreover, for any symplectic basis {e^, hj;j = 1 ,..., s} 


+yjPj) 

i=i 


on T>,where R{ej) = qj, R{hj) = and [a:i, pi,..., are coordinates of 

vector 2 ) in the basis. Here the canonical observables qj,Pj]j = 1,... ,s are 
selfadjoint operators in Ri satisfying the Heisenberg CCR on P 


[qj^Pk] c i5jki, [qj^qk] = o, [Pj,Pk] = o. (26) 

In physics the symplectic space is the phase space of the classical system (such 
as electro-magnetic radiation modes in the cavity), the quantum version of 
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which is described by CCR. Then s is number of degrees of freedom, or 
“normal modes” of the classical system. 

The state given by density operator p in "H is called Gaussian, if its 
quantum characteristic function 

= TipWiz) 

has the form 

((){z) = exp (^im{z)-^a{z,z)^ , (27) 

where m is a real linear form and a is a real bilinear symmetric form on Z. 
A necessary and sufficient condition for fl27p to define a state is nonnegative 
definiteness of the (complex) Hermitian forrr0 a (z, z') — | A (z, z') on Z or, 
briefly: 

a > 1 a. (28) 

We will agree that the matrix of a bilinear form in fixed a symplectic base 
is denoted by the same letter, then (|28il can be understood as inequality for 
Hermitian matrices, where a is real symmetric and A is real skew-symmetric. 

A Gaussian state is pure if and only if a is a minimal solution of this 
inequality, see e.g. [38]. Operator J in Z is called operator of complex 
structure if 

= -I, (29) 

where I is the identity operator in Z, and the bilinear form A(z, Jz') is an 
(Euclidean) inner product in Z, i.e. 

A{z,Jz') = A{z',Jz){=-A{Jz,z')); (30) 

A{z, Jz) >0, z & Z. (31) 

The following characterization can be found in [16], [39] : 

Proposition 5 The minimal solutions of the inequality are in one-to- 
one correspondence with the operators J of complex structure in Z given by 
the relation ^ 

a {z, z') = 

^ A complex-valued real-bilinear form J (z, z') on Z will be called Hermitian if /3 (z', z) = 

dOW). 
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In this way to every complex structure corresponds the family of pure 
Gaussian states fl27p with different values of m which are called the J-coherent 
states. The state with m = 0 is called J- vacuum. Let po be a vacuum, then 
any associated coherent state is of the form W{z')poW{z')*, as follows from 
the relation 

W{z’)*W{z)W{z') = e^p[-i^{z,z')]W{z) 

and from nondegeneracy of the form A{z, z') due to which m{z) = A{z, 
Operator S' in Z is called symplectic if A{Sz, Sz') = A{z, z') for all z, z' G 
Z. The unitary operators W{Sz) satisfy the CCR hence by the Stone- 
von Neumann uniqueness theorem there is a unitary operator Us in "H such 
that 

W{Sz) = U*sW{z)Us, zeZ. 

The map S —)■ f /5 is a projective representation of the group of all symplectic 
transformations in Z, sometimes called “ metaplectic representation” [2] as it 
can be extended to a faithful unitary representation of the metaplectic group 
which is two-fold covering of the symplectic group. 

Similarly, T is antisymplectic if A{Tz,Tz') = —A{z,z') for all z,z' G Z. 
There is an antiunitary operator Ut in Ti such that 

W{Tz) = UtW{z)Ut, zeZ. 


Let Za,Zb be two symplectic spaces with the corresponding Weyl sys¬ 
tems. Consider a channel <1> : A —> B. The channel is called Gaussian if 
the dual channel satisfies 


^*[Wb{z)] = WA{Kz)exp 




z G ZB^ 


(32) 


where K : Zb ^ Za is a linear operator, I a linear form and /i is a real 
symmetric form on Zb- In terms of characteristic functions of states. 


(j)B{z) = (l)A{Kz)exp 



1 

2 


■p{z,z) 


It follows that Gaussian channel maps Gaussian states into Gaussian states. 
A converse statement also holds true |16j . 

A necessary and sufficient condition on parameters {K, I, p) for complete 
positivity of the map $ is (see c® nonnegative definiteness of the Hermitian 
form 

2 ), z' —)■ p {z, z') - ^ [As {z, z') - Aa {Kz, Kz')] 
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on Zb, or, in matrix terms (if some bases are chosen in Za, Zb), 

/i > ^ [As - K^^aK] , (33) 

where * denotes transposition of a matrix. The proof using explicit construc¬ 
tion of the representation of type ([7]) is given in [TT], see also [39]; below in 
Proposition [TTl below we give such a construction for an important particular 
class of Gaussian channels. 

We call the Gaussian channel extrem^ if /i is a minimal solution of the 
inequality fl3^ . This terminology stems from the fact that the minimality of 
fi is necessary and sufficient for the channel $ to be an extreme point in the 
convex set of all channels with hxed input and output spaces [38] . 

Additivity hypothesis for quantum Gaussian channels: The ad¬ 
ditivity properties (lT5il . (1T6|1 hold for any pair of Gaussian channels <hi, $ 2 - 

Hypothesis of quantum Gaussian minimizers: For any function 
f E d the inhmum in fl22|) is attained on a pure Gaussian state p. 

Any Gaussian channel has the covariance property 

^Wa{z)pWa{z)*] = Wb{K^z)^p]Wb{K^z)* (34) 

where is the symplectic adjoint operator dehned by the relation 


Ab {K‘'za, Zb) = Aa {za, Kzb) ■ 

It follows that the value Tr/(d>[p]) is the same for all coherent states W{z)pqW{z)* 
associated with a vacuum state po- 

These two problems turn out to be closely related. In what follows we 
describe positive solution for both of them in a particular and important 
class of Gaussian channels with gauge symmetry. However both conjectures 
remain open for general quantum Gaussian channels. 

3.2 Complex structures and gauge symmetry 

Given an operator of the complex structure J one dehnes in Z the Euclidean 
inner product j{z,z') = A{z,Jz'). Then one can dehne in Z the structure 

® In quantum optics one speaks of quantum-limited channels |20] . 
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of s—dimensional unitary space Z in which iz corresponds to Jz and the 
(Hermitian) inner product^ is 

j(z,z') = + = y^j{z,z') -ij{z, Jz')]. 

From fl2^ . fl5U]) it follows that J is symplectic, that is A{Jz,Jz') = 
A{z, z') for all z,z' ^ Z . With every complex structure one can associate the 
cyclic one-parameter group of symplectic transformations {e‘^'^;(p G [0,27r)} 
which we call the gauge group. Hence, by the Stone-von Neumann uniqueness 
theorem, the gauge group in Z induces the one-parameter unitary group of 
the gauge transformations G [0,27r)} in % according to the formula 

W{e^-^z) = U;W{z)U^. (35) 

For the future use it will be convenient to introduce the complex parametriza- 
tion of the Weyl operators by dehning the displacement operators 

D{z) = W{Jz), z G Z. (36) 

A state p is gauge invariant if p = U^ppU* for all <p, which is equivalent 
to the property (f){z) = z) of the characteristic function. In particular, 
Gaussian state fl27)l is gauge invariant if m[z) = 0 and a {z, z') = a {Jz., Jz'). 
By introducing the Hermitian inner product in Z 

« (z, z') = {z, z') - ia (z, Jz')], 

we have ct (z, z) = (z, z) since a (z, Jz') is skew-symmetric; moreover, the 

condition fl28p is equivalent to nonnegative dehniteness of the Hermitian form 
OL (z, z') — ij(z, z') on Z : 

a > Ij. (37) 

This follows from application of the following Lemma to the form 

z' —^ /3(z, z') = a{z, z') - ^ A ( 2 ;, z'). 

The relation fl37l) can be considered as the inequality for the matrices of the 
form, provided a basis is chosen in Z. In an orthonormal basis, j = I is the 
unit matrix. 

®In accordance with convention accepted in mathematical physics, the inner product 
is complex linear with respect to z' and anti-linear with respect to z. 
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Lemma 6 Let /3{z,z') be a bilinear complex-valued Hermitian form on real 
vector space Z, satisfying f3{Jz, Jz') = (3{z,z'), where J is a linear operator 
such that .P = —I. Then f3{z,z') is nonnegative definite i.e. 

'^CjCk/3{zj,Zk) >0 (38) 

jk 

for any finite collection {zj} C Z and any {cj} C C, if and only if 

Re/5(^, z) ± Im/9(^, Jz) > 0 for all z E Z. (39) 

Proof. fl59D ^> fl55D : We have I3{z,z') = Rel3{z,z') + i\m.l3{z,z'), where 
1111 / 9 ( 2 :, z') is skew-symmetric, hence \t!i(3{z,z) = 0. By using the fact that 
f3{Jz,z') = —l3{z,Jz') we obtain that also Re(3{z, Jz') is skew-symmetric, 
hence Re(3{z, Jz) = 0. Thus 

Ref3{z^ z) ± \m.f3{z^ Jz) = l3{z, z) =F i(d{z, Jz). 

Now introduce complexification z ■(-)■ z by letting Jz ■(-)■ iz and define two 
Hermitian forms on the complexification Z of Z : 

f3'^{z,z') = l3{z,z')Tifd{z,Jz'). (40) 

Then f3~ is sesquilinear i.e. complex linear with respect to z' and anti- 
linear with respect to z, while is anti-sesquilinear. From (ITOD . 1 1391) . 

/3^(z, z) = Re/9(^, z) ± Im/ 9 ( 2 ;, Jz) > 0 for all z G Z, (41) 

hence by (anti-)sesquilinearity 

{zj,Zk) > 0 . 

jk 

By adding the two inequalities corresponding to plus and minus, we get fl55]) . 

Conversely, fl38|) ^^ fl39|) : Applying fl38|) to the collection {zj.,Jzj} C 
Z, {cj, ±icj} C C we obtain 

CjCk [/3 {zj, Zk) ± i(d {zj, Jzk)] > 0, 

jk 

hence the forms fITOl) are nonnegative definite. By (anti-)sesquihnearity of 
these forms, this is equivalent to fITTD i.e. (1591) . ■ 
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Assume that in Za, Zb operators of complex structure Ja, Jb are fixed, 
and let U^, be the corresponding gauge operators in 'Ha,'Hb acting ac¬ 
cording ([35]). Channel $ : A —)■ 5 is called gauge-covariant , if 

4[f/> (U;^)’] = U^i\p] (U^y (42) 

for all input states p and all 0 G [0, 2tt]. For the Gaussian channel (13^ with 
parameters (iC, I, p) this reduces to 

/(^) = 0 , KJb-JaK = 0, p{z,z') = p{Jbz,Jbz'). 

The relation fl3^ for gauge-covariant Gaussian channel takes the form 

<F*[T)s(z)] = T>a(Kz) exp [-^(z, z)], z G Z^, (43) 

where 

M [js - K^JaK] (44) 

The equivalence of fl44|) and (l3^ is obtained by applying the lemma [6] to 
the Hermitian form 

f5{z, z') = p{z, z') - ^ [Ab{z, z') - Aa{Kz, Kz')] . 

Ghannel <F : A —)■ i? is called gauge-contravariant, if 

mp W)’] = (uyy ( 45 ) 

for all input states p and all 0 G [0, 2tt]. For the Gaussian channel fl52]) with 
parameters (A', I, p) this reduces to 

l{z)=Q, KJb + JaK = 0, p{z,z') = p{Jbz, Jbz'). 

The relation fl32p for gauge-contravariant Gaussian channel takes the form 

<F*[T)b(z)] = T)A(-AKz)exp [-/x(z,z)], z e Zb, (46) 

where A is antilinear operator of complex conjugation, A^ = J, A^ = —A in 
Za such that AJa + JaA = 0, and K = —AK is complex linear operator from 
Zb to Za- Here 

/x>±^[js + K*jAK]. (47) 

The last condition is obtained by applying Lemma [6] to the Hermitian form 

f3{z,z') = p{z,z') - ^[Ab{z,z') - Aa{Kz,Kz')] 

= '-[Ab{z,z') + Aa{Kz,Kz')]. 
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3.3 Attenuators and amplifiers 

In what follows we restrict to channels that are gauge-covariant or contravari- 
ant with respect to hxed complex structures. Therefore, to be specihc, we 
consider vectors in Z as s—dimensional complex column vectors, where the 
operator J acts as multiplication by i, the corresponding Hermitian inner 
product is j(z,z') = z*z' and the symplectic form is A{z,z') = 2Imz*z', 
where * denotes Hermitian conjugation. The linear operators in Z commut¬ 
ing with J are represented by complex s x s— matrices. The gauge group acts 
in Z as multiplication by Gaussian gauge-invariant states are described 
by the modihed characteristic function 

0(z) = TrpZi)(z) = exp {—z*olz) , (48) 

where a is a Hermitian correlation matrix satisfying ck > 1/2 as follows from 
([27]). For the given complex structure, the unique minimal solution of the 
last inequality is |l, to which correspond the vacuum state po and the family 
of coherent states {pz; z G Z} , such that pz = T)(z)poT)(z)*. One has 

TrpwT)(z) = exp f2ilmw*z “ 



where |zp = z*z. 

Let Z/i, Zs be the input and output spaces of dimensionalities s^, sb- We 
denote by sa = dimZ^, sb = dimZ^ the numbers of modes of the input 
and output of the channel. The action of a Gaussian gauge-covariant channel 
fl43|l can be described as 

^*[Db{z)] = DAiKz)exp{-z*^z), zeZb, (49) 

where K is complex sb x sa— matrix, /x is Hermitian sb x sb— matrix satis¬ 
fying the condition (see [30] ) 

AX > (1b - K*K ), (50) 

where Ib is the unit sb x ss—matrix. This follows from fl4T|) by taking into 
account that the matrix of the form j(z,z') in an orthonormal basis is just 
the unit matrix I of the corresponding size. Later we will need the following 
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Lemma 7 The map ( [7^ is injective if and only z/0 KK* > 0 (in which 
case necessarily sb > sa)- 

Proof. Injectivity means that <h[pi] = implies pi = p 2 - But <h[pi] = 

$[^ 2 ] is equivalent to Trpi$*[Zi)s(z)] = Trp 2 $*[-DB(z)], i.e. TrpiZi)A(Kz) = 
Trp 2 -DA(Kz) for all z G 'Lb- By irreducibility of the Weyl system, this prop¬ 
erty is equivalent to RanK = Za, i.e. KerK* = {0} or KK* > 0. ■ 

The channel (14^ is extreme if /r is a minimal solution of the inequal¬ 
ity (1501) . Special cases of the maps are provided by the attenuator 
and amplifier channels, characterized by matrix K fulhlling the inequali¬ 
ties, K*K < I and K*K > I respectively. We are particularly interested in 
extreme attenuator which corresponds to 

K*K < Is, /X = ^ {Ib - K*K), (51) 

and extreme amplifier 

K*K >1b, ^ 1 ( K*K - 1b) . (52) 

Denoting by z the column vector obtained by taking the complex con¬ 
jugate of the elements of z, the action of the Gaussian gauge-contravariant 
channel fHbll is described as 

r[DB(z)] =DA(-K^)exp(-zVz), (53) 

where /x is Hermitian matrix satisfying the inequality 

/x>i(lB + K*K), (54) 

which follows from (071). Here z is the column vector consisting of complex 
conjugates of the components of z. These maps are extreme if 

^^ = ^ilB + K*K). (55) 

The following proposition generalizes to many modes the decomposition 
of one-mode channels the usefulness of which was emphasized and exploited 
in the paper [20] (see also [13] on concatenations of one-mode channels): 

^ For Hermitian matrices M, N, the strict inequality M > N means that M — N is 
positive definite. 
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Proposition 8 Any Gaussian gauge-covariant channel ^ : A ^ B is a 
concatenation <l> = <l> 2 ° of extreme attenuator : A ^ B and extreme 
amplifier ^2 '■ B ^ B. 

Any Gaussian gauge-contravariant channel ^ : A ^ B is a concatenation 
of extreme attenuator $1 : A —)■ i? and extreme gauge-contravariant channel 
^2 ■ B —)■ B. 

Proof. The concatenation = $ 2 °*^*! of Gaussian gauge-covariant channels 
$1 and $2 obeys the rule: 


(56) 

(57) 


K = KiKj, 

ft = Kj/JiKj + ftj. 


By inserting relations 



into fl57|) and using fl56|l we obtain 



(58) 


from the inequality fl50l) . By using operator monotonicity of the square root, 
we have 


IK 2 I > Is, IK 2 I > |K 
The hrst inequality (j58ll implies that choosing 



(59) 


and the corresponding P 2 = ^ — Ib) , we obtain extreme ampliher 


^2 ■ B ^ B. 


Then with 



we obtain, taking into account the second inequality in (l58|l and also Lemma 
[9] below. 



which implies K^Ki < I^, hence Ki with the corresponding - i (ifi - K*Ki) 
give the quantum-limited attenuator. 
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Lemma 9 Lei M > K*K, i/ien KM K* < where means (generalized) 
inverse. 


Proof. By the definition of the generalized inverse, 


u = 


\u v\ 


sup 


:'ueRanM,-!;^0 n*Mn 


By inserting in place of u and using Cauchy-Schwarz inequality in the 
nominator of the fraction, we obtain 




sup 

E Ran M, u0 


u*uv*K.*K.v 


u u. 


In the case of contravariant channel the relations (156|1 . (j57jl are replaced 
with 


K = K 1 K 2 , (62) 

M + Ma- (63) 

By substituting 

,.i=1(i-k:k,), ^2 = i (K;K2 + Ib) 

into fl63|) and using fIMl) we obtain 

|K2|' = K;K2 = m + ^(K*K-Ib) >K*K. (64) 

Taking K 2 =|K 2 |, /X 2 = |(|K 2 P + Ib) gives extreme gauge-contravariant 
channel $2 : -B —)■ -B . With 


Ki = K IK 2 


(65) 


we obtain, by using Lemma [H 


KiK* 


K(|K2r)'K* 
1 


K 


fi + -[K*K 



K*<1a 




( 66 ) 


which implies KiK^ < 1^, with the corresponding give the extreme at¬ 
tenuator $1 : A —)■ 5. ■ 
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Remark 10 In the case of gauge-covariant channel, the eguality in m) 
shows that KK* > 0 implies KiK^ > 0, while the ineguality /x (K*K - Ib) 
implies KiK^ < 1a- In the case of gauge-contravariant channel, the inegual¬ 
ity fi > ^ (Ib + K*K) implies 0 < KiK^ < 1a via ^6E) . 

Proposition 11 The extreme attenuator with matrix K and extreme atten¬ 
uator with matrix K = \/Ta'^^KK* are mutually complementary. 

The extreme amplifier with matrix K and gauge-contravariant channel 
with matrix K = \/KK* — 1a are mutually complementary. 


Proof. For the case of one mode see [13] or [33], Sec. 12.6.1. We sketch 
the proof for several modes below. Dehne Zg ~ Z^, — Zb, so that 

Z = Z^ ©Zb — Zb©Zb 

In the case of attenuator consider the block unitary matrix in Z : 


K VIa - KK* 

VIb - K*K -K* 


(67) 


which dehnes unitary dynamics U in TL = TLa ®TLd — TLb ® He by the 
relation U*Dbe{'^be)U = Dad{^'^be)- Here zbe = [zb zb]*, Dbe{2be) = 
Db{ Zb) ® De{ze), and the unitarity follows from the relation 

K^/1b - K*K = (68) 


Let pe) = pq be the vacuum state, pA = p svn arbitrary state. Then the 
formulas (Cj), ([H]) dehne the mutually complementary extreme attenuators as 
described in the hrst statement. The proof is obtained by computing the 
characteristic function of the output states for the channels. For the state of 
the composite system pbe = U{p® Pd)U* we have 


0be(zbe) — Trf/(p © Pd)U* [ZIb(zb) ® De{^e)] 

= Tr(p © pb)H* [T)b(zb) ®-Db(zb)] H 

= Tr(p©pB) Z1 a(Kzb + Kzb) © T)b(\/Ib - K*Kzb - K*z^) 


= 0a(Kzb + Kzb) exp 
By setting ZB = 0orzB = 0we obtain 


--|v/Ib-K*Kzb-K*zb| 


(69) 


0i?(zB) = 0A(KzB)exp 

<(>e(zb) = 0A(KzB)exp 


--z^ (Ib-K*K)zb 


--zJ^KK*z^ 
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as required. 

In the case of amplifier, set 


V 


K -VKK* - IaA 

-AVK*K - Is AK*A 


where A is the operator of complex conjugation, anticommuting with mul¬ 
tiplication by i. By using the property A{Az,Az') = —A{z,z'), we obtain 
that V corresponds to a symplectic transformation in Z generating unitary 
dynamics U in Ti. Let again po be the vacuum state of the environment. 
Then the formulas ([7]), ([8]) dehne the mutually complementary channels as 
described in the second statement of the Proposition, and the proof is similar. 

To show that V is a symplectic transformation, introduce the matrices 


1b 0 

, Vi = 

K 

VKK* - Ia 

1- 

0 

1 

> 

[ VK*K - 1b 

K* 


1b 0 
0 — 


Notice that V = 0Vi©, and Vj'SVi = H, which means that Vi preserves 
the indehnite Hermitian form (t{zbej ^'be) ~ ^E^'e- By taking into 

account 


Abe{®zbei ®z'be) — ^b{zb, z'q) — Ae{ze, z'e) — lm.a{zBEi z'be)i 


we obtain 


ABE^y Zbe-,^ z'ee) — 1®ZbEi^ 1®z'ee) — 1"^(z{@ZbEi®z'ee) — ^Be{zbEi z'be)i 


as required. ■ 

Again, later we will need the following 

Lemma 12 Let \ A ^ B he an extreme attenuator with 0 < KiKj; < I^, 
then ‘hi[Tp] = Pp> (a pure state) if and only if Pp is a coherent state. 

Proof. According to Proposition [TTl the complementary channel $1 is an 
extreme attenuator with the matrix K = such that 0 < K < 

I^. Its output is also pure, d’liPp] = as the outputs of complementary 
channels have identical nonzero spectra by Lemma [H Thus 
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where ^ is the vacuum vector and U is the unitary operator in "H 
implementing the symplectic transformation corresponding to the unitary 
fl67|l in Za®Zd ~ Zb®Ze, with Zd — Zb, Ze ^ Za Denoting by 


0(z) = TrP^DA(z), 0 '(zb) = TrP^/Ds(zs), (j)E{zE) = TrP^^ DEi^E) 

the quantum characteristic functions and using the relation ([69]), we have 
the functional equation 


0'(zb)0£;(ze) = 0(KiZb+Kze) exp 


-IiVIb-KJKiZb - K\ie 


. (70) 


By letting z^; = 0, respectively z = 0, we obtain 


0'(zs) 

0_e(ze) 


0(KiZb) exp 
0 (Kz£;) exp 


--|VIb-K*KiZb| 


1 , 


-^iK^z^r 


thus, after the change of variables z = KiZg, z' = Kz^, and using ([68]), the 
equation (1701) reduces to 

0(z)0(z') = 0(z + z') exp [Rez*z']. 

The condition of the Lemma ensures that RanKi =RanK = Za- Substi¬ 
tuting a;(z) = 0(z) exp [||zp] , this becomes 

ci;(z)ci;(z') = a;(z -I- z') (71) 


for all z,z' G Za- The function a;(z), as well as the characteristic function 
0(z), is continuous and satisfies uj{—7.) = a;(z). The only solution of flTTl) 
satisfying these conditions is the exponent C(;(z) = exp [Rmw*z] for some 
complex w. Thus 


0(z) = exp 


Rm w*z 


1 

2 


z 


2 


is the characteristic function of the coherent state Pw/ 2 - ■ 
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3.4 Gaussian optimizers 

The following basic result for one mode was obtained in |53] . Here we present 
a complete proof in the multimode case, a sketch of which was given in |22]. 


Theorem 13 (i) Let ^ be a gauge covariant or contravariant channel and 
let f be a real concave function on [0,1], such that /(O) = 0, then 

Tr/($[p]) > Tr/(<h[p^]) = Tr/(<l>[po]) (72) 


for all states p and any coherent state (the value on the right is the 
same for all coherent states by the unitary covariance property of a Gaussian 
channel [34\))- 

(a) Let f be strictly concave, then equality in |7l| ) is attained only if p is 
a coherent state in the following cases: 

a) sb = sa and $ is an extreme amplifier with pi = \ (K*K — 1b) > 0; 

b) sb > Sa, the channel is gauge-covariant with KK* > 0 and 


pi>- (K*K - 1b] 


(73) 


c) Sb > Sa, the channel is gauge-contravariant with KK* > 0 and 
A. > i (Is + K*K). 


Proof, (i) We hrst prove the inequality fl72l) for strictly concave /. Then the 
inequality for arbitrary concave / follows by the monotone approximation 
f{x) = lim^io/£(a^), since fe^x) = f{x) — ex‘^ are strictly concave. Also, by 
concavity, it is sufficient to prove fl72|) only for p = P^j,. 

By Proposition [HI $ = $ 2 ° where $1 : A —)■ H is an extreme at¬ 
tenuator and $2 : H H is either extreme amplifier or extreme gauge- 
contravariant channel. Any extreme attenuator maps vacuum state into vac¬ 
uum. Indeed, 


[po] L)b{z) 


Trpo*hi[T)s(z)] 

Trpo7^A(Kz) exp (Is - K*K) z 

exp = TtpoDb{z). 
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Therefore Tr/(<h[po]) = Tr/($ 2 [po])- Then it is sufficient to prove fl72D for all 
extreme amplifiers and all extreme gauge-contravariant channels $ 2 . Indeed, 


assume that we have proved 

Tr/($2[Pp]) > Tr/(<h2N). (74) 

for any state vector tlj. Consider the spectral decomposition <hi [P^] = PjPt /,^, 
where pj >0, then 

W('S>1P*1) = TV/(42|4i|P*]]) (75) 

> EP2Tr/(1’2[P«J) (76) 

j 

> ^(42^01) (77) 

= lV/(42|4i|p„||)=TV/(4b„|). (78) 


Then, according to the second statement of Proposition [TT] and Lemma [2] 

W(42|7’*1) = TV/(i2l-P*]). 

where $2 is an extreme amplifier and $2 is an extreme gauge-contravariant 
channel. Thus it is sufficient to prove fl74|) only for an extreme amplifier 
$2 : 5 —)■ 5, with Hermitian matrix K 2 > Is- 

The following result is based on a key observation by Giovannetti. 

Lemma 14 For an extreme amplifier $2 : S —)■ 5, with matrix K 2 > Is, 
there is an extreme attenuator such that for all tf 

4>2(Pp)~(4>2 0$;)(Pp). (79) 

Proof. By Proposition ITT] and Lemma |2] <h 2 (Pp) ~ l> 2 (Pi/i) for all fj E Pb, 
where <i >2 is extreme contravariant channel with the matrix K = a/K^ — Is- 
Dehne the transposition map T : P —)■ P by the relation T[P(z)] = 
P(—z). The concatenation $ = To $2 is a covariant Gaussian channel: 

4>*[P(z)] = !>; or[P(z)] = P(^Ki-lBz)exp . 

Applying decomposition from Proposition [HI namely the relation ( l58l) . gives 
$ = $ 20 $'^, where $2 is the original ampliher, and : P —)■ P is another 
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extreme attenuator with matrix Ki = This implies the relation 

(uni)- ■ 

Lemma [m and Lemma 0] imply 

Tr/($ 2 (Pp)) = Tr/ (($2 o $;) (Pp)) . (80) 

Again, consider the spectral decomposition of the density operator 

p'>0. 


By concavity, 


Tif ((4^ o 4'i) [P*]) > ^P'W (42|P*,]) 


( 81 ) 


Since / is assumed strictly concave, then p —)■ Tr/(<h 2 [p]) is strictly con¬ 
cave m- Assuming that Pp is a minimizer for the functional flSOp . we con¬ 
clude that ^ 2 [Ppj] must all coincide, otherwise the above inequality would be 
strict, contradicting the assumption. From Lemma[7]it follows that Pp. = Ppi 
for all j and for some ip' ^Pb-, hence, assuming that Pp is a minimizer, the 
output *hi[Pp] = Pp> is a pure state. 

Since Ki = the condition of Lemma IT^ is fulhlled if K 2 > Is. 

In this case, if Pp is a minimizer, the Lemma implies that Pp is a coherent 
state. Thus we obtain the inequality fTMl) for the ampliher $2 with K 2 > Is 
and strictly concave / . In this way we also obtain the case a) of the “only 
if” statement (ii). 

In the case of ampliher $2 with K 2 > Is, we can take any sequence 
> 1b, ^ K 2 , and the corresponding amplihers $ 2 "^- Then 

Tr/($ 2 "'^[p]) —)■ Tr/(<F 2 [p]) for any concave polygonal function / on [0,1], 
such that /(O) = 0, and any p G ©("Ha)- This follows from the fact that any 


h”’ip]-42ip] 


-)■ 


such function is Lipschitz, \ f{x) — f{y)\ < M:\x — y\, and 

0. It follows that holds for all extreme amplihers $2 in the case of 
polygonal concave functions / . For arbitrary concave / on [0,1] there is 
a monotonously nondecreasing sequence of concave polygonal functions fm 
converging to / pointwise. Passing to the limit m —)■ cx) gives the inequal¬ 
ity fl7^ for arbitrary extreme ampliher, and hence, fl72|) holds for arbitrary 
Gaussian gauge-covariant or contravariant channels. 
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(ii) The “only if” statement in the cases b), c) are obtained from the 
decomposition $ = $ 20 $^ and the relations d7S]) - d77D by applying argument 
similar to the case of extreme ampliher. Notice that the conditions on the 
channel $ imply that in the decomposition $ = $2 o $1 the attenuator 
$1 is dehned by the matrix Ki such that 0 < KiK^ < Ia (see Remark 
[TU]) . Applying the argument involving the relations flSUll - flST]) with strictly 
concave / to the relations fl75]) - fl75]l . we obtain that for any pure minimizer 
Fp of Tr/(<h[F^]) the output of the extremal attenuator ‘hi[Fp] is necessarily 
a pure state. Applying Lemma fT8l to the attenuator $1 we conclude that Fp 
is necessarily a coherent state. ■ 

3.5 Explicit formulas and additivity 

Proposition 15 For any p > 1 and any Gaussian gauge-covariant or con- 
travariant channel 



= (Tr$[po]")'/", 

(82) 

RvW 

= Fp(<h[po]), 

(83) 

am 

= Hi^po]), 

(84) 


where po is the vacuum state. 

The multiplicativity property a holds for any two Gaussian gauge- 
covariant (contravariant) channels $1 and $ 2 , as well as the additivity of 
the minimal Renyi entropy l[T5\) and of the minimal von Neumann entropy 

m- 

Proof. The hrst statement follows from Theorem [13] by taking f{x) = —x^, 
/(x) = —xlogx. 

If <hi and $2 are both gauge-covariant (contravariant), then their tensor 
product $1 (g) $2 shares this property. The second statement then follows 
from the expressions (152]) - (IHTI ) and the product property of the vacuum 
state Po = ^ ® p^q \ which follows from the dehnition. ■ 

From the dehnitions of gauge-co/contravariant channels (H9l) . ( l53l) . it 
follows that the state ^[po] is gauge-invariant Gaussian with the correlation 
matrix /x -|- K*K/2. The spectrum of <F[po] is computed explicitly leading to 
the expressions [H] 

= [det [(p + K*K/2 + 1b/2Y - (p + K*K/2 - 1^/2)"]]-'/" 
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and 


^($) = tr^(/x+ (K*K - I^) /2), (85) 

where g{x) = (a: + 1) log(x + 1) — xlogx and tr denotes trace of operators 
in Z. In the last case we nsed the formnla for the entropy of Ganssian state 

dH]) HO]: 

H{p) = tig{cx - 1/2). 

We now tnrn to the classical capacity of the channel <I>. In inhnite dimen¬ 
sions, there are two novel featnres as compared to the sitnation described in 
Sec. 12.41 First, one has to extend the notion of ensemble to embrace continnal 
families of states. We call generalized ensemble an arbitrary Borel probabil¬ 
ity measnre tt on The average state of the generalized ensemble vr is 

dehned as the barycenter of the probability measnre 

Ptt = j p'K{dp). 

6 (Ha) 

The conventional ensembles correspond to hnitely snpported measnres. 

Second, one has to consider the inpnt constraints to avoid inhnite valnes 
of the capacities. Let F be a positive selfadjoint operator in Ha, which 
usnally represents energy in the system A. We consider the inpnt states with 
constrained energy: TrpF < F, where F is a hxed positive constant. Since 
the operator F is nsnally nnbonnded, care shonld be taken in dehning the 
trace; we pnt TrpF = Xdmp{X), where mp{X) = TrpF(A), and F(A) is 
the spectral fnnction of the selfadjoint operator F. Then the constrained y— 
capacity is given by the following generalization of the expression fITTI ): 

Cx{^,F,E)= snp y(7r), (86) 

7v\TrpT^F<E 

where 

y(7r) = F(4>[p^]) - j H{<l>[p])n{dp) (87) 

6(Ha) 

To ensnre that this expression is dehned correctly, certain additional condi¬ 
tions npon the channel $ and the constraint operator F shonld be imposed 
(see |3n], Sec. 11.5), which however are always fnlhlled in the Ganssian case 
we consider below. 
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Denote = F(8)/---0/ + -- - + J(8)---(8)/(8)F, then the constrained 
classical capacity is given by the expression 

C'($,F,E)= (88) 

n^oo n ^ 

Now let $ be a Gaussian gauge-covariant channel, and consider gauge- 
invariant oscillator energy operator F = where e = [e^fc] is 

a Hermitian positive dehnite matrix, Oj = -^{qj + ipj) ~ the annihilation 
operator for j-th mode. For any state p satisfying TrpF < oo, the hrst 
moments Trpoj and the second moments Trpa*afc, Trare well dehned. 
For gauge-invariant state Trpoj = 0 and Tr^aja^ = 0. For a Gaussian gauge- 
invariant state fHHjl 

a-1/2= 

see e.g. [33] . 

Proposition 16 The constrained classical capacity of the Gaussian gauge- 
covariant channel $ is 

C{^;F,E) = C^{^;F,E) 

= max trc/(KVK + /x+(K*K-I b)/ 2)-trc/(/x+(K*K 

u: tTue<E 

The optimal ensemble tt which attains the supremum in consists of co¬ 
herent states pz = Da{z)pqDa{z)*, z G distributed with gauge-invariant 
Gaussian probability distribution Q^{d'^^z) on Z^ having zero mean and the 
correlation matrix u which solves the maximization problem in [8^) . 

Proof. Gonsider a Gaussian ensemble tTi, consisting of coherent states pz = 
D^(z)po-Da(z)*, z E Za, with gauge-invariant Gaussian probability distribu¬ 
tion Qjj{df^z) on Za having zero mean and some correlation matrix u. It is 
dehned by the classical characteristic function 

J exp (2ilmw*z) Q^{d‘^’^w) = exp {—z*uz). 

■Za 

By using the covariance property flM|) of Gaussian channel, we have 
Hi^pz]) = if(<F[D^(z)poD^(z)*]) = ih(<F[po]) = trp(/x + (K*K - I^) /2), 


(89) 

Ib)/2). 
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which does not depend on z, and hence it gives the value of the integral term 
in fl87p . Integration of the characteristic functions of coherent states gives 

Trp^^DA(z) = exp (-z* {u + 1^/2) z) . 

Then u = and Trp^^F = = true. 

The state *h[p 7 r^] is gauge-invariant Gaussian with the correlation matrix 
K* {jy + Ia/ 2) K+^^, hence it has the entropy tr 5 f(K*^'K+^t+(K*K — I^) /2). 
Thus for the Gaussian ensemble 

X(7r,) = tr ^(KVK + /x + (K*K - 1^) /2) - + (K*K - I^) /2). (90) 

Summarizing, we need to show 

C{^;F,E) = C^{^;F,E)= sup xM- (91) 

u: tri>e<E 

Let us denote by Q the set of Gaussian gauge-invariant states in I-La- 


Lemma 17 


max iL (<h®"') < n max H(^\p]). (92) 

pW;TrpWFW<n£ ^ Lr J / 


Proof. We hrst prove that 


sup i/(<h®"[p(’")]) < n sup i7(<h[p]). 

p(”)-.Tr p(")F(^)<nE p:TrpF<E 


(93) 


Indeed, denoting by pj the partial state of in the j—th tensor factor 

of "H®” and letting p = ^ have 


//(4®“|p(")|) < < nH{m), 

1 = 1 

where in the hrst inequality we used subadditivity of the quantum entropy, 
while in the second - its concavity. Moreover, TrpF = iTrpGlFl"-) < E, 
hence fl95]) follows. 

Using gauge covariance of the channel <h, we can then reduce maximiza¬ 
tion in the right hand side of fl93|) to gauge-invariant states. Indeed, for a 
given state p, satisfying the constraint TrpF < E the averaging 

1 

Pav = ^ U^pU*d(p 
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also satisfies the constraint, while H{^[p]) < H{^[pav]) by concavity of the 
entropy. 

Finally, we nse the maximum entropy principle which says that among 
states with hxed second moments the Gaussian state has maximal entropy 
(see e.g. [39], Lemma 12.25). This proves fl9^ . ■ 

We have 


max H(^\p]) 

p:peg,TrpF<E 


max tr^(KVK + /x+(K*K-I b)/2)]. (94) 

u:tr:ue<E 


Now let v be the solution of the maximization problem in the righthand side. 
To prove (1591) observe that 


nxM < nG^(d>,F,E) < nE) 

< max iL (<|)®"- — min 

pW:Trp(")F(")<n£; ^ ^ p(") 


if (<1.®^ [pW]) , 


By using Lemma [T7| and Proposition [15] we see that this is less than or 
equal to 


n 


max F ($ [p])-H (<F [po]) 

p:pGy,TrpF<E 


= nxM, 


where the equality follows from fl9T|) and fl90|) . 

Thus C'^(<1>®”, FG)^ 77,F) = nCy^{^, F, E) and hence the constrained clas¬ 
sical capacity (1551) of the Gaussian gauge-covariant channel is given by the 
expression (159]) . ■ 

Similar argument applies to Gaussian gauge-contravariant channel (153|) . 
giving the expression (159|) with e replaced by e. Indeed, in this case the state 
<F[pjr^] is gauge-invariant Gaussian with the characteristic function 


Tr<I)[p^^]F(z) = exp (-(Kz)* {u + 1^/2) Kz - z*puz) 
= exp (—(Kz)* (P-I-Ia/2) Kz — z*iiz), 


with the correlation matrix K* (P -i 1^/2) K -i /x. On the other hand, trize = 
trPe, so that redehning as u, we get the statement. 

The maximization in (l89l) is a hnite-dimensional optimization problem 
which is a quantum analog of “water-hlling” problem in classical information 
theory, see e.g. usmo]. It can be solved explicitly only in some special cases, 
e.g. when K, n, e commute, and it is a subject of separate study. 
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3.6 The case of quantum-classical Gaussian channel 

Consider affine map which transforms qnantnm states p G &{'H) into prob¬ 
ability densities on Z 


p pp(z) = TrpD(z)po^(z)*, (95) 

where -D(z) are the displacement operators, po is the vacnnm state with the 
qnantnm characteristic fnnction 


0o(z) = Trpo^(z) = exp 

The fnnction Pp(z) is bonnded by 1 and is indeed a continnons probability 
density, the normalization follows from the resolntion of the identity 

/ D(z)po/^(z)*^ = 1 . 

Proposition 18 Let f be a concave function on [0,1], such that /(O) = 0, 
then for arbitrary state p 

r p2s r Ms 

/ > / /(Ppw(z))^. (96) 

Jz TT Jz 71 

Proof. For any c > 0 consider the channel <Fc defined by the relation 

= J Tr[pT>(c-^z)po/^*(c"^z)] p^. (97) 

The map fl^ is a Gaussian gauge-covariant channel such that 
<F*[Zl(z)] = D{cz) exp 


(c^ + 1), 



cf. |23] . Therefore by Theorem [T3l 

'IV/(4,1p|) > T1/(4,K]) (98) 

for all states p and any coherent state p^. We will prove the Proposition ITSl 
by taking the limit c —)■ cxd. 
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In the proof we also use a simple generalization of the Berezin-Lieb in¬ 
equalities |H]: 


/(P(z))^ < Tr/(a) < f f{p{z)) 

TT J7. 




TT* 


(99) 


valid for any quantum state admitting the representation 

a= p(z)pz 




71 ° 


with a probability density p(z). In the right side of p(z) = Trcrpz. In 
the inequalities (19^ one has to assume that / is dehned on [0, cx)) (in fact, 
p(z) can be unbounded). We shall assume this for a while. 

Taking a = <hc[p], from fl97)) we have 


p{z) = 4- TTpD{c ^z)poD*{c ^z) = -^Ppic ^z) 


while 


p(z) = Trp^^c[p] = / p{w)Trp^p, 




7T° 


( 100 ) 


We use the well-known formula, see e.g. 

TrpzPw = exp[-|z - w^]. 

By introducing the probability density of a normal distribution 


^2s 


^c(z) = 

tending to <5—function when c —>■ 00 and substituting this into fllUUp . we have 

p(z) = J p(w) gi(z — w) 

= [ d^^w' Pp(w') gi(z - cw') 


= -^Pp*(lc{c z). 


(101) 


With the change of the integration variable c ^z ^ z, the inequalities 
become 


/(c ‘^"Pp{z))^<c 2"Tr/(<I)Jp]) < [ /(c ^^Pp*qc{z)) 


d""z 


71° 
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Substituting p = p^, we have 


[ /(c <c ^*Tr/(<h,[pw]) < [ f{c qdz))^. 

Jz " Jz " 

Combining the last two displayed formulas with (IMll we obtain 

j giPpiz))^- J gipp^iz))^ 

> [ £/(Pp(z))^- [ g{pp*qc{2))^, (102) 

Jz ^ Jz ^ 

where we denoted g{x) = f{c~‘^''^x), which is again a concave function. More¬ 
over, arbitrary concave polygonal function g on [0,1], satisfying 5 f( 0 ) = 0, 
can be obtained in this way by dehning 

= j a;e[0,c-2'*] 

\ g{l) + g'{l){x — xe[c“^^,oo) ’ 

hence fll02p holds for any such function. Then the right hand side of the 
inequality fll02p tends to zero as c —)■ oo. Indeed, for polygonal function 
\g{x) — g{y)\ < >c\x — y\, and the asserted convergence follows from the 
convergence Pp * t Pp in Li : if p(z) is a bounded continuous probability 
density, then 

lim / \p * qc{z) — p{z)\(f^z = Q. 

C-^QO 

Thus we obtain (lOOP for the concave polygonal functions /. But for arbi¬ 
trary continuous concave / on [0,1] there is a monotonously nondecreas¬ 
ing sequence of concave polygonal functions /„ converging to / . Applying 
Beppo-Levy’s theorem, we obtain the statement. ■ 

4 Appendix 

Consider a gauge-covariant channel <I> such that the matrices K*K and 
commute (in particular, this condition is satished by extreme amplihers and 
attenuators). These channels are diagonalizable in the following sense. We 
have 

K = V^K.Vb, m = , 


40 












where V^, are unitaries and K^, are diagonal (rectangular) matrices 
with nonnegative values on the diagonal. Then K*K = V^K^V^, and 

$[p] = UbMUapUX]U*b, (103) 

where Ua, Ub are canonical unitary (“metaplectic” 0) transformations act¬ 
ing on Ha, Hb such that 

U*^Db{z)Ub = Db{Vbz), U\Da{z)Ua = Da{Vaz), 

To describe the action of “ diagonal” channel in more detail, we have to 
consider separately the cases sa = sb, sa < sb and > sb- 
In the case = ss we have 


Krf = diag = diag • 

Then where, in self-explanatory notations, 

= DjikjZj) exp {-fij \zjf) . 

In the case sa < sb 


K,= 


diag 

0 


(104) 


where 0 denotes block of zeroes of the size (sb — sa) x Then 

where for j = 1,..., the one-mode channels are given by ( I104p . and 
^[sA+iv,ss] -g vacuum state of the modes + 1, • • •, sb- 
In the case sa > sb 

Kd = [ diag 0 ] 

where 0 denotes block of zeroes of the size sb x (s^ — sb), and 

^d[p] = [Tr,^+i,...,,^p], 

where Tr^g+i^...^^^ denotes partial trace over the last sa — sb modes of the 
operator p. 

There is a similar reduction to the diagonal form for gauge-contravariant 
channels. 
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